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Abstract 


The temperature observed (~ 10^ °K) for the hydrogen atoms of 
interstellar HI clouds is shown to be a consistent measure of the 
Coulomb interactions of cosmic rays. These interactions could also 
give rise to a high electron concentration (~ 10"^ cnT^). The expected 
intensity of Hey emission by the hydrogen atoms excited by cosmic rays 
is found to be comparable to but less than that measured from the 
galactic disk. The inner bremsstrahlung X-radiation by electrons 
ejected from interstellar hydrogen atoms is indicated as a direct 
measure of the spectrum of the cosmic rays responsible for the ionization. 
A detailed study of the non-stellar X-rays from Hi clouds is recommended 
for the determination of the spatial and spectral distribution of the 
low energy (<* 10® e .v. /nucleon ) cosmic rays that propagate throughout 
the galactic disk. 
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Introduction 

Cosmic ray measurements during the recent period of minimum 
solar modulation indicate that there is an appreciable flux of low 
energy highly ionizing cosmic ray particles in the interstellar 
medium (see Fig. 1). However, the modulation (Balasubrahmanyan , 

Boldt and Palmeira 1967) appears even at solar minimum to be suffi- 
ciently strong and complicated to make a precise inference of the 
interstellar spectrum quite unlikely unless there is the input of 
some additional relevant information directly from outside the 
solar system. This probelm is here approached from the outside by 
considering the Coulomb interactions of these particles in the inter- 
stellar medium and examining the possible electromagnetic signatures 
that could be measured from observations made from within the solar 
system. 

It is important to note that most of the interstellar gas is 

in HI clouds of neutral hydrogen (for a review see Dieter and Goss, 

1966). The Doppler broadening of the 21 cm line of atomic hydrogen 

indicates a gas temperature of — 10^°K for such clouds. The possible 

relevance of low energy cosmic rays for maintaining this temperature 

was noted by Hayakawa (1960). Recent work by Field (1962) and 

Smith (1966) on the cooling mechanisms for these clouds indicates 

/ 1 

that the characteristic cooling time (< 4 x 10 years ) is 
much shorter than the meantime (~ 7 x 10 years ) between the 
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cloud-cloud collisions described by Kahn (1955) as the events 
responsible for most of the heating. Hence, a steady-state input 
such as the heating by cosmic rays now appears compelling, and 
the problem is here re-examined within the context of new information. 


Thermal Balance 

The energy dissipated by cosmic ray nuclei (E > E Q ) in the 
interstellar medium via Coulomb interactions with the constituent 
hydrogen atoms is given per gram-second by: 


q <E 0 ) 


4tt 


f fdE\ /dJ\ dE 

J \dxj \dE ) T\ 


e.v. 

g-sec. 


( 1 ) 


where: E = kinetic energy per nucleon, 

T| = fraction of cosmic ray nucleons that appear as 
free protons, 

dE ^ e.v. loss/ (g-cm‘2 ) 5 for protons traversing atomic hydrogen, 
dx 


dJ = interstellar protons . 
dE (cm^-sec-sr-e.v. ) 

The Coulomb interaction of cosmic rays with interstellar 
hydrogen atoms results in 

i) the ejection of atomic electrons, 

ii) inner bremsstrahlung x-rays generated by the ejecting 
electrons, and 

iii) the excitation of atomic levels. 
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The ejected electrons slow down by heating the gas, and finally 
recombine with ions to produce characteristic hydrogen emission 
lines. The heated hydrogen atoms cool by collisional excitation of 
infra-red emitting levels of other atomic species, especially 
oxygen (Smith, 1966). 


Under steady-state conditions, the conservation of energy 
appears as: ^2) 

( \ /Inner \ / Hydrogen \ / \ / > 

Cosmic Ray ) = ( Bremsst rah lung) + ( Line ) + ( Atomic ) + ( Ionic 

Heating J \ x-Radiation J \ Emission / ^Coolingy/ \ C ° oling , 


q(o) 


'IB 


+ 


^ + x 


3/2 (NkT ) 4- 3/2 (NkT ) 


where q^ = e.v. / (g-sec. ) that appears as inner bremsstrahlung 


X-radiation > 


= e.v./ (g-sec. ) that eventually emerges as the spectral 
line (hv^) , 

Tj = relaxation time (sec.) for cooling by the (j)th atomic species, 

t+ — relaxation time (sec.) for cooling by collisional 

excitation of several ionic species via intermediary 
electrons , 

kx = 10 -2 e.v. (for T« 10 2o K), 

23 

N = Avagadro Number (6 x 10 ) . 
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The X-radiation from inner bremss trahlung is the only out- 
put of the cloud that directly tracks the spectrum of the cosmic 
ray nucleons, but as far as the energy balance (?) is concerned 
the X-rays make an entirely negligible contribution. As described 
by Hayakawa, Nishimura, and Takayanagi (1961), the important hydrogen 
atomic emission lines in this situation are Ly-a, lb, and lb. 

o 

Smith (1966) has estimated that near T ^ 10 °K, the major contri- 

bution to cooling by atomic collisional excitation is accounted 
for by the interstellar abundance of oxygen. The relaxation time 
(t c ) for this process is not a sensitive function of temperature in 
the vicinity of 10 °K and is given by 

T 0 = 12 x 10^- sec -» 4 x 10^ years. (3) 

- 1 -3 

For the case where the electron concentration is high (n £ > 10 cm ) , 
cooling by the Seaton (1955) mechanism of ionic collisional excitation 
by electrons becomes important. To provide for this we use a value 
of t_|_ extracted from the work of Takayanagi and Nishimura (1960), viz 

t+ = 10^ sec. ( 4 ) 

n e 

2 

The ion-electron density (n e ) produced in cool (~10 °K) 

interstellar Hi regions (Hayakawa, 1960) by cosmic rays is fixed by 
q(o)rf For a steady-state situation in a medium of pure hydrogen 
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the production via the ionization process is balanced by recombination 


as follows: 


CT r V e n e = 


where: a r = 2 x lCr^ cm^ = recombination cross section (Bethe and 

Salpeter, 1957) at T = 10^ °K, 

V e = thermal electron velocity at 10^ °K, 
n^ = Hi hydrogen number density (~ 10 cm ), 

AE = 36 e.v. — mean energy dissipated by a cosmic ray per 
ion pair formed (Rossi, 1952), 

N = Avagadro Number 
From (5) we obtain at 10^ °K that 


-7 1/2 _3 

7 x 10 7 q(o ) cm J 


The insertion of this expression (6) for into (4) gives 
directly in terms of q(o), viz 


= 10 J 


q (o ) "1/2 seCi 


From the work of Hayakawa (1963) and Hayakawa, Nishimura and 
Takayanagi (1961), we estimate 


q TR = 3 x 10° q ( 0 ) 


( 8 ) 
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Y, q i = 4 X 10- 1 q ( 0 ) 


(9) 


The substitution of (3), (7), (8), and (9) into (2) 
yields an equation for q(o) alone, with the solution 


q (o ) - 3 x 10 ^ e.v. / (g-sec) . 


( 10 ) 


This thermal solution (10) for q(o) may be used to obtain 
the corresponding thermal solution for n e , via (6), as 


10 


-1 


cm 


-3 


(ID 


This value for n e is likely to be an upper limit since there 
might be other recombination processes competitive to simple 
electron-proton attachment. On the basis of the low frequency 
(< 5 x 10^/sec) fall-off for the galactic radio spectrum in directions 
towards the galactic poles, Hoyle and Ellis (1963) conclude that 
<">>• the average electron density in our region of the galactic 
disk, is given by 


<^n e ^> » 1CT 4 T 3 / 4 . 


( 12 ) 


If the electrons responsible for the low frequency galactic radio 
absorption by the galactic disk reside mainly in Hi clouds then 
T ^ 10 °K for these electrons. That Hi clouds occupy ~ 10% of 
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the volume of the galactic disk would imply that n 

e 

whereby (12) gives 



n e « 10 <n e ^> » 10" 1 cm" 2 3 * * 6 . 


(13) 


Radiation 

The electromagnetic radiation from HI clouds is a measurable 
consequence of the cosmic ray phenomenon described here. The 
intensity (1^) of various spectral lines is given by 
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where (S) is the optical thickness in hydrogen atoms/cm for 
penetrating radiation (e.g. H>); for diffuse radiation (e.g. Ly-<y) 

(S) is comparable to the thickness of the galactic disk, in 
hydrogen atoms/cm . 

Using the results of Hayakawa, Nishimura and Takayanagi (1961) 
for q^/(E q^) with the approximation (9), q(o) given by (10), S = lO^/cm^ 
for Ly-a, and S = 10^2 f 0 r H , we obtain 

h 2 

10 e.v./(cm <*sec-sr) . (15) 

and 

6 2 

Ijk, =“ 10 e.v./ (cm -sec-sr) , 


( 16 ) 
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This value (16) for galactic Hey radiation corresponds to ~2 
Rayleighs, which is within an order of magnitude of the Hey intensity 
associated with the galactic disk, (~ 5 -* 20 Rayleighs) as measured 
by Dufay, Dufay and Huu-Doan ( 1 96 1 ) • 


From the work of Hayakawa (1963) we may estimate that, for inner 
bremsstrahlung , 


(17) 


dl ~ 3 x 10‘ 9 (S)_ 

dhv 4 tt 



500 (hv) 


e . v. 


(cm^ sec-sr-e.v.) 


Substitution of the definition (1) for q into (17) yields an 
integral transform for the interstellar cosmic ray proton spectrum 
(dJ/dE), viz 



(3 x 1 0~ V ) S 


oc> 

/(§)(§) 

(500 h v ) 


(18) 


dE 

r 


(cn? 


e.v. 

-sec-sr-e.v. 


) 


Since the ionization loss rate (dE/dx) is known as a function of E, 
a determination of (dl/dhv) for the inner bremsstrahlung X-rays 
gives a direct representation of the cosmic ray proton spectrum. 

Consider the X-ray spectral intensity at hv & 30 keV ; this 
corresponds to Eo = 15 MeV/nucleon for the interstellar cosmic ray 
spectrum. As described in the next section (see Table 1), a 
consistent interstellar spectrum gives q(Eo) 3 x 10 10 e.v./(g-sec) 
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for Eo — 15 MeV/nucleon. The spectral intensity of 30 keV 

X-rays for directions that traverse most of the galactic disk is 

23 ? 

obtained by inserting S « 10 atoms/cm z into (17) and is 
evaluated as 


( ~ 1 9 , , e - v - 

\ dhv/ (cm z -sec-sr-e. v. ) 

30 keV 


(19) 


This intensity (19) is comparable to the observed background 
(Bleeker et al., 1966) of X-rays at ~ 30 keV. It is important to 
note that a marked anisotropy is expected for the X-rays discussed 
here since the intensity viewed normal to the galactic plane would 
be diminished by an order of magnitude. Measurements of the non- 
stellar X-rays from HI clouds should yield the spatial distribution 
and spectrum of the low energy (< 10^ MeV/nucleon) cosmic rays that 
propagate throughout the galactic disk. 

Interstellar Cosmic Ray Spectrum 

Recent measurements (Balasubrahmanyan et al., 1966; Comstock 
et al., 1966; Fan et al., 1966; Biswas et al., 1966; Meyer et al., 
1967; Hormann and Winckler, 1966; Ormes and Webber, 1966) of the 
spectrum for cosmic ray , H^, , He^ , He^, C^, and 0^ 


nuclei indicate that : 
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i) the interstellar spectrum for cosmic ray alpha part- 
icles is similar in form to that for protons, on a 
kinetic energy per nucleon basis (Balasubrahmanyan et al., 
1967), over the band of maximum observed flux (i.e. ^ 

10® -» 10^ e.v. /nucleon) , 

ii) the mean cosmic ray traversal through the interstellar 
medium is ~ 6 g/cm , but could be as low as 3 g/cm at 
the upper and lower energy extremes of the spectrum 
(Fan et al., 1966; Biswas et al., 1966; Meyer et al., 
1967), 

iii) the spectral shape for C f N, and 0 nuclei (Comstock 

et al., 1966) relative to alpha particles at low energies 
(< 500 MeV/nucleon) may be explained (Cowsik et al., 1966) 
by a steady-state model whereby the interstellar path 
distribution is exponential, favoring short paths, and 

iv) the source spectrum for cosmic ray protons and alpha 
particles cannot flatten at low energies, as would be 
the case for a power law in total energy (Meyer et al., 
1967). 


A proton source spectrum that is compatible with present observations 
is : 



(E/m 0 c ) 


-2.5 


m 0 c 


(unit E-cm -sr-sec ) 


( 20 ) 
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where E is the kinetic energy and m 0 c is the rest energy. 

This source spectrum (20) is to be modified by the energy 
degradation of the protons that traverse an exponential distribution 
(Cowsik et al., 1966) of interstellar paths (X), the calculation 
has been made for both ^X^ ~ 6 g/cm^ and <^X^ = 3 g/cm^. These 

spectra are exhibited in Figure 1, along with a curve given by 
Gloeckler and Jokipii (1967) as a summary of proton measurements 
made during 1965, a period of minimum solar modulation. The ratio 
of the interstellar spectrum (dJ/dE)^ for X^> = 6 g/cm^ 

to the spectrum observed in 1965 is plotted as a function of 1/g in 
Figure 2. Over the band 50 -> 800 MeV/nucleon, where solar modulation 
is significant, we note that the ratio may be approximated by 


(dJ/dE ) x 

(dJ/dE ) 1965 


exp (2.7/g). 


( 21 ) 


This formula (21) for the solar modulation of cosmic rays over the 
indicated energy band is of the form shown by Balasubrahmanyan , 
Boldt and Palmeira (1967) to be an adequate description during 
several years near the period of minimum solar activity. 


The interstellar cosmic ray proton spectrum discussed here 
( (dJ/dE ) x ), exhibited in Figure 1, may be used to estimate a 
value for q consistent with cosmic ray particle observations. 
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The contributions to q from the energy losses by other nuclei 
(mainly He^) is here approximated (Balasubrahmanyan et al., 1967) 
by taking r\ ~ 1/2 in equation (1). Using the known energy degrad- 
ation of protons traversing atomic hydrogen for energies above E 0 
(15 MeV) we obtain from (1) that q for the interstellar spectrum 
proposed here (see Fig. 1) is: 

q(E 0 ) = (2 - 4) x 10 10 e.v./g-sec. (22) 

for ^x^ ~ (3 -6) g/cm^ 

Shifting E 0 by 10 MeV/nucleon changes this value of q by less than 
507 o . As shown in Table 1, the interstellar spectrum discussed by 
Gloeckler and Jokipii (1967) gives a comparable value of q for E q 
= 15 MeV. However, the latter spectrum steepens sharply below 50 Mev/ 
nucleon and this has the effect of changing q by over an order of 
magnitude when E 0 is lowered by < 10 MeV . 

The value of q obtained (22) on the basis of this spectral 
analysis of observed cosmic rays matches the value (10) for the 
thermal solution. Therefore, it appears consistent to attribute 
the heating of Hi regions mainly to the Coulomb interactions of 
cosmic rays. Models of the interstellar spectrum are compared for 
this role in Table 1; a previous extreme estimate based on a 
minimal solar modulation (Balasubrahmanyan et al., 1966; Brunstein 
and Cline, 1966) is seen to be inadequate (see Table 1, last entry). 
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A measure of the electromagnetic energy deposition by cosmic 
rays during a lifetime is given by 

q x\^ — q <x> (e.v./cnf*). (23) 

As exhibited in Table 1, this is comparable to the kinetic energy 
density for the interstellar proton spectrum examined here (Fig. 1). 
Most of this kinetic energy density is accounted for by those protons 
of high ionizing power; the kinetic energy density for the relativistic 
protons (E >500 MeV ) of this interstellar spectrum is less than 

o 

1 e.v./cm . The power U dissipated by all cosmic rays via Coulomb 
interactions with the entire mass (Ginzburg and Syrovatskii, 1964) 
of galactic gas (~ 3 x 10^ 2 g) ± s given by: 

U = Mq = (3 x lO^ 2 g) (2 x 10^ e.v./(g-sec) ) (24) 

= 6 x 10 52 e.v./sec. 

-» 4 x 10^® ergs/sec. 

This value of U for the dissipated power is equal to the total power 
discussed by Ginzburg and Syrovatskii (1964) as required for the 
steady-state replenishment of cosmic rays that escape from the galaxy. 
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F igure Captions 


Figure 1. The Observed" curve is a summary of the proton spectrum 
as given by Gloeckler and Jokipii (1967) , based upon 
balloon and satellite borne observations by several lab- 
oratories (Balasubrahmanyan et al., 1966; Comstock et al., 
1966; Fan et al., 1966; Biswas et al., 1966; Meyer et al., 
1967; Hofman and Winckler, 1966; Ormes and Webber, 1966) 
during 1965. The curves ”<[3 g/cm2^>" and M <6 g/cm^^> ,! 
represent the indicated power law source spectrum after 
traversing interstellar matter via an exponential dis- 
tribution of paths with a mean of 3 g/cm^ and 6 g/cm^ 
respectively. 

Figure 2. The dotted curve is the ratio of the interstellar spectrum 
for <6 g/cm^]> as shown in Fig. 1 to the spectrum observed 
in 1965 (Fig. 1). The solid line is exp(2,7/g). g is 
the proton velocity in units of the velocity of light (c) . 
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